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I. INTRODUCTION 

Research in plasmonics is expanding its domains into 
several sub-fields due to significant advances in ex- 
perimental techniques .'i'^ The unique optical proper- 
ties of the surface plasmon-polariton (SPP) resonance, 
being the very foundation of plasmonics, find intrigu- 
ing applications in optics of nano-materials,'^lt2l materi- 
als with effective negative index of refraction,'i2Hl2! jj. 
rect visualization,^'^ photovoltaics,'^^'^ single molecule 
manipulation,'iSH23 and biotechnology.'2iM231 Theoretical 
modeling of optical properties of metal nanostruc- 
tures is conventionally based on numerical integra- 
tion of Maxwell's equationSjI^^ltiSl although simulations 
within time-dependent density functional theory ap- 
peared recently for small atomic clusters .I^SEi] More- 
over, current theoretical models are quickly advanc- 
ing toward self-consistent simulations of hybrid materi- 
als: metal/semiconductor nanostructures optically cou- 
pled to ensembles of quantum emitters.f^^l This method- 
ology, based on numerical integration of corresponding 
Maxwell-Bloch equations, brings new insights into nano- 
optics as it allows for the capture of collective effects. 

The molecular optical response in a close prox- 
imity of plasmonic materials is greatly enhanced by 
SPP modes leading to the discovery of the single 
molecule spectroscopy.l^StSl] Recently, experiments per- 
formed on current carrying molecular junctions started 
to appear.!^^!!^ Theoretical modeling of molecule-SPP 
systems utilizes the tools of quantum mechanics for the 
molecular part. In particular, studies of optical response 
of isolated molecules absorbed on metallic na nopar ti- 
cles utilize Maxwell-Bloch (Maxwell-SchrodingerjPliUfUl 

equations or near fie ld-tim e dependent density functional 
theory formulations .I^SSni 

Realistic molecular devices are open quantum systems 
exchanging energy and electrons with surrounding envi- 
ronment (baths). This is especially important in stud- 
ies of molecules in current carrying junctions interacting 
with external fields.®' Usually in such studies the electro- 
magneticfEM) field is assumed to be an external driving 
force! ^^"^^ ^ Recently we utilized the nonequilibrium Green 
function technique to study the transport and optical re- 



sponse of a molecular junction subjected to external EM 
field taking into account near-fields driven by SPP local 
modes, specific for a particular junction geometry.l^^Eni 
Although the formulation allows us to describe the molec- 
ular junction with formation of the local field by SPP 
excitations in the contacts taken into account explicitly, 
the molecular influence on formation of the local EM field 
was disregarded in these studies. Note that such influ- 
ence was shown to have measurable effects in plasmonic 
spectrum.i22lllIllIlS! 

When a molecule located near metal surface is driven 
by a strong EM field, one can expect to observe signifi- 
cant changes in the total EM field due to radiation emit- 
ted by the molecule. Such radiation although quickly 
degrading with the distance from molecular position can 
nevertheless noticeably alter the local EM field. Since the 
latter is driving the molecule, transport characteristics of 
the junction may be significantly modified. This calls for 
a self-consistent treatment, where both SPP excitations 
and molecular response participate in formation of the 
local EM field. 

Here we extend our previous considerations by tak- 
ing into account complete electrodynamics and molecu- 
lar junction response in a self-consistent manner com- 
bining Maxwell's equations with electron transport dy- 
namics. The molecule is treated as a pointwise source 
in the Ampere law. We demonstrate the importance 
of the molecular response in the formation of the local 
field for an open molecular system far from equilibrium. 
The effect is shown to be important for proper descrip- 
tion of the junction transport characteristics. The paper 
is organized as follows. Section |lT] presents a transport 
model of the molecular junction. Section III describes 



the methodology of computing the EM field taking into 
account molecular response. The results are presented in 
section HVl Section IVl summarizes our work. 



II. 



MOLECULAR JUNCTION SUBJECTED TO 
EXTERNAL EM FIELD 



We consider a junction with a molecular bridge {M) 
connecting between two contacts (L and R). The bridge 
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FIG. 1: (Color online) A sketch of the junction. 



is formed by D two-level systems with the levels repre- 
senting ground [g) and excited [x) states of the molecule. 
Each of the two level systems is subjected to a classical 
local EM field E{t) (see section 



III 



for details of its cal- 
culation). Electron transfer is allowed along the chain of 
ground (excited) levels of the bridge. The contacts are 
taken in the form of bowtie antennas, and are assumed to 
be reservoirs of free electrons each in its own equilibrium 
with electrochemical potentials /i^ and /x^j, respectively 
(see Fig. [ij. The Hamiltonian of the system reads (here 
and below e = h = 1) 



H{t)=HM{t)+ J2 (Hk + Vk 



K=LM 



(1) 



s=g,x 
D 



D~l 



m— 1 



drriKS + H.C 



keK s=g,x 
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where HM{t) and Hk are Hamiltonians of the molecu- 
lar bridge (M) and the contacts [K — L,R), and Vk is 
coupling between them. In Eqs. ([2])-(j4| {dms) and 
c|, (cfe) are creation (annihilation) operators for an elec- 
tron on the level s of the molecular bridge site m and 
state k of the contact, respectively. Emit) is the local 
time-dependent field at bridge site m, and flms.ms' — 
{ms\fl\ms') is the matrix element of the transition molec- 
ular (vector) dipole operator between states \ms) and 
|ms'). For simplicity below we assume that the transi- 
tion dipole moment is the same for all bridge sites and has 
only one non-zero component, fimg.mx = fJ-gx for any to. 
ts (s — g, x) and Vk^mKS are matrix elements for electron 
transfer in the molecular bridge and between molecule 
and contacts, respectively, and niK = 1 (D) ioi K — L 
(R). Note that treating the external field classically al- 
lows us to account for arbitrary time dependence exactly 
(i.e. beyond perturbation theory).^ 
We follow the formulation of Ref. |311 Time-dependent 



current at interface K (L or R) is 



I Kit) 



-IniTr 



(^G<(t,t) + I ^Me) G^{t,e)^ 



(5) 

where Tr[. . .] is a trace over the molecular subspace, 
fxie) = [e^""'"^)/'^ + is the Fermi-Dirac distribu- 
tion in contact K, is the molecular dissipation matrix 
due to coupling to contact K 

(e) EE 27r E K^lSl,fcVfc,,„2S2^(e - e^), (6) 



miSi, 7712^2 



and G^*^""^ is a matrix in the molecular basis of the lesser 
(retarded) projection of the single particle Green func- 
tion, defined on the Keldysh contour a^^ 



G 



(n, T2) = -i{Tc dm.si (n) dLs2 ('^2)) (7) 



Here Tc is the contour ordering operator and ri^2 are the 
contour variables. In Eq.([5]) G''(t, e) is the right Fourier 
transform of the retarded projection of the Green func- 
tion ^ 



G'-(i,e) = J dt' e'<'-''^ G^t,t') 



(8) 



Note that in Eq.([5]) and below we assume the wide band 
limii|£5l in the metallic contacts. 

The Green functions in ^ satisfy the following set of 
equations of motion^SHH 

z|G'-(t, e) = I - (^el - HM(t) + G^t, e) (9) 

tj^G<{t,t) = [UM{t);G<{t,t)] - '-{T;G<{t,t)} 



E / ^fK{e){r'<G^{e,t)-G^{t,e)r^) 



(10) 



K=L,R 



where I is the unity matrix, Hj\,/ (t) is a representation of 
the operator ([2| in the molecular basis, T = ^k=l r ^ 
[...;...] and^...;...} are the commutator and anti- 
commutator, and G"(e,t) = [G''(i, e)]^. The first or- 
der differential equations ^ and ( [lO| are solved starting 
from the initial condition of the biased junction steady- 
state in the absence of the optical pulse, E{t = 0) = 



Gl{e)^G'-{t = 0,e) = 
G< EEG<(t = 0,i = 0) 



eI-HM(t = 0) + -r 



- E 

K=L,R 



de_ 
2^ 



G^„{e)r'<fK{e)G^„{e) 



(11) 
(12) 



where Gg(e) = [G5(e)]^ 

Below we calculate the charge pumped through the 
junction by the optical pulse 



Qit)= f dt 
Jo 



,lL{t')~lR{t') 



-ht 



(13) 
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where lL,R{t) are defined in Eq.([5|, and Jq is the steady- 
state current 

lo^l ^Tr[r^G5(6)r«Gg(6)](/L(6)-/^(6)) (14) 



III. SELF-CONSISTENT ELECTRODYNAMICS 

The time evolution of electric, E, and magnetic, H, 
fields is considered according to the set of Maxwell's 
equations (written here in SI units) 



dH(r,t) ^ , 



(15a) 



eo^^^ - V X i?(f, t) - J(r, i), (15b) 

where /io a-nd eo are the magnetic permiability and di- 
electric permittivity of the free space, respectively, and 
J{t) is the electric current density. Note that magneti- 
zation is disregarded in Eqs. (15a I and (15b), since we 



assume both molecule and contacts to be non-magnetic. 

A molecule located at site m ( = fv„) and driven 
by local electric field E{f„i,t), yields time-dependent re- 
sponse, which enters Ampere's law as a polarization cur- 
rent density 



J{rm,t) 



dPmjt) 

dt 



(16) 



where S is the Dirac delta-function. The polarization 
depends on molecular characteristics through the molec- 
ular density matrix, which in turn is affected by the local 
field. In our model two-level systems of the molecular 
bridge (2]) are assumed to occupy sites of the FDTD grid. 
Molecules contribute to the polarization at their site ac- 
cording to 



Pmit) = 2Im 



mx.mg '^mg.m 



(17) 



The resulting system of coupled differential equations, 
Eqs (15l-(15b), is solved simultaneously with EOMs for 
the Green functions of the quantum system, Eqs. 
(10). The Maxwell's equations are discretized in time 



and space and propagated using the finite-difference 
time-domain approach (FDTD).'^ We employ three- 
dimensional FDTD calculations utilizing home-build par- 
allel FORTRAN-MPI codes on a local multi-processor 
cluster.ISSl In spatial regions occupied by a plasmonic 
nanostructure (a bowtie antenna in our case) we employ 
the auxiliary differential equation method to account for 
materials dispersion. The dielectric response of the metal 
is modeled using a standard Drud e form ulation with the 
set of parameters describing silver .I^SEH The Green func- 
tions EOMs are propagated with the fourth order Runge- 
Kutta scheme. 

Within described self-consistent model the local elec- 
tric field Emit) = E{rm,t) in Eq.(|2| driving a molecular 
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FIG. 2: (Color online) Map of the instantaneous electric field 
strength, [El{r,t) + E'y{r,t) + El{r,t)Y^'^, at a distance of 10 
nm from the molecule (the plane is parallel to xy) calculated 
(a) without and (b) with the molecular response. The dis- 
tribution is shown for t — 77.8 fs and 81.7 fs for (a) and (b) 
respectively. See text for parameters. 



junction is thus defined by both SPP excitations in the 
contacts and the local molecular response. In the next 
section we show that the molecular contribution changes 
junction transport characteristics drastically, and in gen- 
eral can not be ignored. 



IV. NUMERICAL RESULTS 

Here we present result of numerical simulations demon- 
strating the importance of a self-consistent treatment of 
the local EM field dynamics. Previous studies consid- 
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ered th e influen ce of an isolated molecule on plasmon 
transfer fSE^^niolecular features in absorption- ^^^" ^-^lini 
and RamarPsMIlll spectra of molecules attached to 
nanoparticles. Below we discuss how molecular junctions 
and electron transport are influenced by a local EM field 
and vice versa in a self-consistent manner. 

Unless otherwise specified parameters of the calcula- 
tions are T = 300 K, Ex = —Eg = 1 eV, = tg = 
0.05 eV, figx = 32 D, Tf^ .^ = T^. = 0.1 eV and 



^12^,12; = ^Dg,Dg — ^-Ol cV (othcr elements of the dis- 
sipation matrix are zero). These parameters are chosen 
to represent a molecular junction with a strong charge- 
transfer trans itiorP ^, and are similar to our previous 
considerations.lSSllll The Fermi energy is taken at the 
origin, Ep = 0, and the bias is applied symmetrically, 
Mi = -^J■R = Vsdf^- 

Following Ref. '48', the incoming incident field is taken 
in the form of a chirped pulse 



So exp 



- i^2)^2 



iioot 



(18) 



where Sq is the incident peak amplitude, luq is the in- 
cident frequency, and (5^ = 2tq/{tq + 4$"^(wo)) and 
fL = — 4$"(cjo)/(tq -I- 4<i>"^(a;o)) are parameters describ- 
ing the incident chirped pulse (tq is the characteristic 
time related to the pulse duration). In the calculations 
below we use £0 = 10^ V/m, = 2 eV, tq = 11 fs, and 
$"(ljo) = 3000 fs2. 

Figure [2] shows instantaneous electric field strength dis- 
tributions in a plane shifted by z = 10 nm parallel to 
xy plane. The distribution is calculated for a junction 
formed by bowtie antennas with single molecule (D = 1) 
placed in the center of the gap. Here Ex — £g = 1.75 eV, 
^ig.lg = ^D.,D. = 0-01 eV, T^x,!. = ^^^3,03 = 0-001 cV, 
and Vsd = 0. Fig. [2^ presents simulations without molec- 
ular response. Fig. [2|d shows results of a calculation 
where both SPP excitations in the contacts and molec- 
ular response are taken into account. One can clearly 
see that even a single molecule drastically changes local 
electric field distribution. 

Sensitivity of the pulse temporal behavior to the molec- 
ular response is presented in Figure [3^. Here a local field 
affected by only SPP modes (dotted line) is compared 
to pulses calculated when molecular response is taken 
into account. The latter may result in both enhance- 
ment (dashed line) or quenching (solid line) of the local 
field depending on the ratio of the pulse frequency, wq, 
to the molecular excitation energy. Ex — Eg. In particular, 
quenching is observed for the laser frequency being below 
the threshold (wq < Ex — £g = 2.25 eV), while frequency 
above the threshold (wq > £x—£g = 1-75 eV) leads to en- 
hancement of the field. To understand this behavior we 
perform a simple analysis treating coupling to the driv- 
ing field as a perturbation, and neglecting the chirped 
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FIG. 3: (Color online) Local EM field at the molecular po- 
sition, (a) Pulse calculated without (dotted line, black) and 
with {ex — eg > LO(i - solid line, red; Ex — Eg < u)o - dashed line, 
blue) molecular response, (b) Maximum local field during the 
pulse vs. molecular excitation energy calculated without (tri- 
angles, black) and with (circles, red) molecular response. See 
text for parameters. 



character of the pulse. This leads to (see Appendix [A|) 



Pi{i) w -£ocos(a;ot) \ljLgx\^ 
Im[G<^ig(e)] . 

Im [Gix,ix{^)\ 



de 
2^ 



(19) 



{Ex ~ Wo) 



e-(e,-cjo)]2 + [ri,aj2]2 

e - (gg + wq) 

[6-(e,-Hwo)]^ + [rig,i3/2]2 



where G*- is the lesser projection of the Green func- 
tion ([t]). Taking into account that in the absence of the 
chirp £'inc(^) — £ocos(wo^) the first term in the right 
side of Eq.([l9| suggests that for populated ground state, 
^ic(ig(^) ~ 1' the molecular polarization oscillates in 
phase with the field for wq < Ex — Eg, and in anti-phase 



for ujq > Ex — Eg. Thus according to Eqs. (15b) and 



(171 the molecular response quenches the field in the for- 
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FIG. 5: (Color online) Current at the left interface as a 



function of time. Shown are differences, A7i = Z^^"^' — , 
between results calculated with I^'"^'', and without, 
molecular response. The calculations are performed for Vsd = 
1.5 V (dashed line, blue) and 2 V (solid line, red). Inset 
shows corresponding difference in charge pumped through the 
junction. See text for parameters. 
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FIG. 4: (Color online) Charge pumped through the junction, 
(a) Difference, AQ = Q'-""'' - Q("°"'=', between results calcu- 
lated with, Q^'"^\ and without, Q^"°'"^\ molecular response 
vs. time for Ex — Eg > uuo (solid line, red) and Ex — £g < ijJo 
(dotted line, blue), (b) Total charge pumped during the pulse 
vs. molecular excitation energy calculated without (triangles, 
black) and with (circles, red) molecular response. See text for 
parameters. 



mer case, and enhances it in the latter. Fig. |3j) illus- 
trates this finding within the exact calculation showing 
the maximum of the total field for different molecular ex- 
citation energies (circles) compared to the maximum of 
the EM field obtained without molecular response (tri- 
angles). Note that contribution of the second term in 



the right side of Eq.(19) is exactly the opposite that of 



the first term, however since the calculations presented in 
Fig. [3] are performed at zero bias, the molecular excited 
state is initially empty, Gf^ ^^{t = 0) ss 0. 

While the local EM field cannot be measured directly, 
it is related to junction characteristics (in particular, its 
transport properties) detectable in experiments. Fig. 
demonstrates the difference in the temporal buildup of 
the charge pumped through the junction, when the 
molecule is considered to be driven by the field obtained 



within the self-consistent model vs. model with only SPP 
excitations taken into account. Initial dip in the charge 
buildup (see dotted line) is related to a time delay of the 
molecule induced pulse for — Eg < ujq (compare solid 
and dashed lines to the dotted line in Fig. |3^). The de- 
lay is caused by the chirped nature of the incoming pulse, 
with initial pulse frequency being lower than the molecu- 
lar excitation energy, which results in suppression of the 
local field at the start of the pulse. Eventually however 
the incoming frequency becomes higher than the molecu- 
lar transition energy. Corresponding enhancement of the 
local field leads to increase in the charge pumped through 
the junction. Note that for — Sg > loq no delay is ob- 
served, and the local field is quenched throughout the 
pulse. Correspondingly effectiveness of the charge pump 
is lower in this case (see solid line in Fig. [4^). 

Figure shows the total charge pumped through the 
junction during the pulse at different molecular excita- 
tion energies. Clearly, the most effective EM field ob- 
tained without the molecular response taken into account 
corresponds to the resonance situation, ujq = Sx — Sg = 
2 eV. When molecular response is included in the model 
the situation is less straightforward. Since local field en- 
hancement is expected for low molecular excitation ener- 
gies, iiJo > Sx — Sg (see Fig. [sj^), the peak in the pumped 
charge distribution is shifted to the left. Note that the 
lower height of the shifted peak is related to the fact that, 
for a lower molecular gap, part of optical scattering chan- 
nels is blocked due to partial population of the broadened 
excited and ground states of the molecule (see Ref.[49lfor 
detailed discussion). 

Note that the importance of molecular response de- 
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FIG. 6: (Color online) Asymmetry in the charge trans- 
fer between positively and negatively chirped incoming laser 
pulses, (5($") — Q(— $"), normalized by their average, Qavg = 
((3(<1?") + (5(— $"))/2. Shown are results calculated without 
(triangles, black) and with (circles, red) the molecular re- 
sponse. See text for parameters. 



pends also on bias across the junction. Indeed, since 
high bias, Vsd > £x —£g, niay inject holes into the molec- 
ular ground state and electrons into the excited state, 
and since populating these states has opposite conse- 
quences for the local field enhancement (see Eq. ( 19 ) 



and the discussion following it), it is natural to expect 
that the molecular response is more important at low bi- 
ases, Vsd < £x — £g- Figure [5] illustrates this conclusion 
with results of our calculations within the self-consistent 
model. Here Tf^ = Ff^ = 0.05 eV. We observe that 
both difference in optically induced current and charge 
pumped through the junction (see inset) is almost negli- 
gible at high biases. Similar reasoning indicates that the 
molecular response at strong incoming fields will be less 
important also due to population of the excited molecular 
state induced by external pulse. 

Asymmetry in the charge pumping relative to the sign 
of the chirp rate was discussed in our recent publica- 
tion (see Fig. 4 in Ref. |48|) . One of the reasons for the 
asymmetry is related to the time spent by the local pulse 
in the region of frequencies at and just below the res- 
onance. This region provides the main contribution to 
charge transfer (see discussion of Fig. 3 in ref. l48|) . Since 
time spent in this region by the positively chirped pulse is 
smaller than that by the pulse with equal negative chirp 
rate (the positively chirped local pulse is shorter), one 
expects to observe an asymmetry as represented by the 
result of calculations using local EM field influenced only 
by SPP modes driving the junction (see curve with tri- 
angles in Fig. [6|. However as discussed above, it is this 
pre-resonance region where molecular response quenches 
local field, thus diminishing (or even overturning) the 
asymmetry relative to the chirp rate sign (see curve with 
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FIG. 7: (Color online) Effect of the self-consistent treatment 
on local field and level population in a 3-sites molecular bridge 
(D = 3) as functions of time. Shown are (a)-(c) local field 
calculated without (dotted line, black) and with (sx — Sg > ijJo 
- solid line, red) molecular response for the three molecular 
sites. Panel (d) shows the difference in population of the 



(nose) 



— /tjg (solid line, blue) and excited, 
(nose) ^^Q-j-^gjj jjjjg^ j.g(^^ states for the first 



ground, Anig 

Anix = n[l^^ - ,.1^ 
molecular site (m = 1). Panel (e) shows the scaled plot of the 
field on the central site (m — 2) for a longer period of time. 
The charge pumped through the 3-sites molecular bridge vs. 
time is shown in panel (f). See text for parameters. 



circles in Fig. |6]). 

Finally, we consider a 3-sites molecular bridge {D = 3) 
to model the spacial nonlocality of molecular polariza- 
tion. Calculations are done for Smx — £mg = 2.25 eV, 
ujQ — 2 eV, and Vsd — 0. Panels (a)-(c) of Fig. [t] com- 
pare the pure plasmonic local field to the field calculated 
when the molecular response is taken into account for 
the three sites of the bridge. Molecular polarization de- 
creases the local field amplitude on the first site, (a), 
and enhances it on the rightmost site, (c). The field at 
the middle site, (b), does not change. The effect can 
be understood following the discussion similar to that of 
Fig. |3] We find that for Sx — £g > ojq increase in pop- 
ulation in the ground (decrease in the excited) levels of 
the molecular sites quenches the local field. Change in 
the populations of the leftmost site, panel (a), resulting 
from self-consistent treatment is shown in Fig. [7[i. We 
see that these changes are in agreement with the corre- 
sponding change in the local field. Similar considerations 
are also hold for panels (b) and (c) (corresponding level 
population are not shown). 

Sclf-consistently calculated electric field on a site in the 
bridge shows a visible beat at large timescale(see Fig. [7^). 
This behavior is related to the Rabi frequency due to the 
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intersite coupling, t^. 

Finally, Fig [7f shows charge transferred through the 
3-site junction as function of time. Decrease in the effec- 
tiveness of the pump is related to quenching of the local 
field on the first site of the bridge, where strong coupling 
to the left contact yields quick resupply of the ground 
level population. Decreased efficiency in pumping the 
charge between ground and excited levels at this site is 
the reason for the overall change in the effectiveness of 
the pump. 



V. CONCLUSION 

We consider a simple model of a molecular junction 
driven by external chirped laser pulses. The molecule 
is represented by a bridge of D two-level systems. The 
contacts geometry are taken in the form of a bowtie an- 
tenna. The FDTD technique is used to calculate the 
local field in the junction resulting from SPP excitations 
in the contacts. Simultaneously we solve time-dependent 
nonequilibrum Green functions equations of motion to 
take into account the molecular contribution to the local 
field formation. 

Note that many works on driven transport assume pure 
incident field to be a driving force acting on the molecule. 
In our recent publications*^ we considered effects of 
local field formation due to SPP excitations in the con- 
tacts on junction characteristics under external optical 
pumping. Here we make one more step by taking into 
account also the molecular response in the driving local 
field dynamics. Within a reasonable range of parame- 
ters we demonstrate that the latter is crucial for proper 
description of the junction transport. We compare our 
results with previously published predictions, and show 
that the molecular contribution may lead to measurable 
differences (both quantitative and qualitative) in charac- 
teristics of junctions. This contribution is especially im- 
portant at low biases and relatively weak external fields 
in the presence of a strong molecular transition dipole. 
In particular, we show that for laser frequencies shorter 
(higher) than the molecular excitation energy the local 
SPP field is usually quenched (enhanced) by molecular 
response. 

Extension of the approach to realistic ab initio calcula- 
tions, taking into account time-dependent bias, and for- 
mulating a methodology for calculations in the language 
of molecular states are the goals for future research. 
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Appendix A: Derivation of Eq.(19) 



To understand trends observed in the exact calcula- 
tions based on Eqs. (9|-(10) and (15a|-(15b), here we 
employ a simple consideration and derive an approximate 



expression for the molecular polarization, Eq.(17), given 



in Eq.(19|. For simplicity we assume that only one pro- 



jection of the molecular dipole is non-zero, and consider 
a single molecule bridge {D = 1). Then the molecular 
polarization is 



(Al) 



Assuming the dissipation matrix, Eq.Q, is diagonal the 
lesser projection of the Green function in Eq.( Al ) is given 
by the Keldysh equation of the form 

Gfg,iAt,t) = Yl f '^^i r dt^Glg^.^tM) (A2) 

s=g.x-'-°° 

X ^fs,ls(^l - ^2) Glg,lsit2,t) 

where 

K=L.R'' 

(A3) 

is the lesser self-energy due to coupling to the contacts. 
We start by neglecting a chirp of the incoming field 



Eincit) = £'ocos(tJoi) 



(A4) 



and treat interaction between molecule and incoming 
field 



^ss' (t^ — ^s' ,s f-^ss ^ijicij^') 



(A5) 



within the first order of perturbation theory. Here s in- 
dicates state opposite to s, i.e. for s = g s = x. 

Within the approximations the retarded Green func- 
tion in Eq. ( A2 ) can be expressed as (similar expression 



can be written for the advanced projection) 



is'(i,^') ~ Ss,s'G[li^{t — t') 



(A6) 



E 



+CX3 



Here G^^-''' is the retarded projection of the Green func- 
tions ([7]) in the absence of external field 

G^iJ'uit - t') = - t')e-^(^-'Ti.,./2)(t-t') 

and d{. . .) is the Heaviside step function. 
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Utilizing (|A4l)-(|A7l) in (|AT])-(|A3| leads to 

Im 

- Im 



Pl{t) « \^J.gx\ 



de_ 
2^ 



--•(OX 



- ujq)] cos{ivot) - [ria;,i^/2] smjujot) 

[e-(£,-c^oF + [ri,.ir./2]' 
- (£g + a;o)] cosjujot) + [rig^ig/2] sin{ujot) 

[e-(£g+C.0)]2 + [rig,lg/2]2 



(A8) 



where we have used the Keldysh equation for the steady 
state situation 



Gls!ls(^) 



Eif=L,fl«/^(e)rfs,i 



(A9) 



[e-e,]2 + [ri,,i,/2]2 
Assuming that detuning is much bigger than levels broad 



enings, \ujo - {e^ - £q)| > Ti^^is (s = g,x), the term 



with sin(u;ot) in (A8) can be ignored. Finally, dress- 
ing the Green functions in Eq. ( |A8| , i.e. taking into 
account diagrams related to population redistribution in 
the molecule due to presence of the driving field, leads to 
Eq.|T9l). 
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